Molecular dynamics study of the thermopower of Ag, Au, and Pt nanocontacts 



O 
<N 

I— » 






B 

o 
o 



m 

o 



^ 



F. Pauly,^ J. K. Viljas,^ M. Biirkle/ M. Dreher,^ P. Nielaba,^ and J. C. Cuevas"^ 

^ Institut fiir Theoretische Festkorperphysik and DFG Center for Functional Nano structures, 

Karlsruhe Institute of Technology, 76131 Karlsruhe, German'^ 

^Low Temperature Laboratory, Aalto University, P.O. Box 15100, 00076 Aalto, Finland 

^ Fachbereich Physik, Universitdt Konstanz, 78457 Konstanz, Germany 

^ Departamento de Fisica Tedrica de la Materia Condensada, 

Universidad Autonoma de Madrid, 28049 Madrid, Spain 

Using molecular dynamics simulations of many junction stretching processes we analyze the ther- 
mopower of silver (Ag), gold (Au), and platinum (Pt) atomic contacts. In all cases we observe 
that the thermopower vanishes on average within the standard deviation and that its fluctuations 
increase for decreasing minimum cross-section of the junctions. However, we find a suppression of 
the fluctuations of the thermopower for the s-valent metals Ag and Au, when the conductance orig- 
inates from a single, perfectly transmitting channel. Essential features of the experimental results 
for Au, Ag, and copper (Cu) of Ludoph and van Ruitenbeek [Phys. Rev. B 59, 12290 (1999)], as yet 
unaddressed by atomistic studies, can hence be explained by considering the atomic and electronic 
structure at the disordered narrowest constriction of the contacts. For the multivalent metal Pt our 
calculations predict the fluctuations of the thermopower to be larger by one order of magnitude as 
compared to Ag and Au, and suppressions of the fluctuations as a function of the conductance are 
absent. 

PACS numbers: 73.63.Rt, 72.15.Jf, 73.23.Ad, 72.10.Fk 



I. INTRODUCTION 

The development of thermoelectric devices for the effi- 
cient conversion of heat into electrical energy or for refrig- 
eration would be an important step towards a more en- 
vironmentally friendly use of energy. Engineered nanos- 
tructures are promising materials in this respect P^ and 
molecular junctions are presently moving into the focus of 
research.^ ^ They seem to be good candidates for achiev- 
ing high thermoelectric figures of merit as a result of 
the discrete energy level structure of the molecule^^ 
and because molecular properties can be controlled by 
chemical synthesis. Already the measurement of the 
thermopower alone provides important insights into the 
electronic structure of molecular junctions, not accessible 
by conductance measurements.^ Thus, electron and hole 
conduction can be distinguished^^ and the alignment of 
molecular levels with respect to the Fermi energy can be 
deter mined. "^"^^J^^ It turns out that the experiments are 
described by a combination of electronic structure and 
transport calculations j ^^ * ^^ * For instance, the linear in- 
crease of the thermopower with molecule length in the 
typical off-resonan t tra nsport sit uation h as been found 
in both experimenlP^ and theory.^^^E^EU 

While molecular junctions appear promising from a 
device-oriented point of view, metallic electr odes are 
typically used to contact the molecules .'^'^^^11 Metallic 
atomic contacts hence constitute important reference sys- 
tems. Until now their thermopower has been studied 
experimentally in the single-atom contact regime only 
by Ludoph and van Ruitenbeek.^ The authors discuss 
mostly their results for Au contacts, but measure ments 
for Ag and Cu contacts have been performed as well.'^^EH 

Motivated by pioneering experiments with high- 



mobility two-dimensional electron gases ,^^^1^ theoreti- 
cal studies have described the thermopower of quan- 
tum point contacts in two and three dimensions us- 
ing the free-electron-gas approximation for adiabatic 
constrict ions .1^2^^ Indeed, such studies have successfully 
predicted maxima of the thermopower at the transition 
between quantized conductance plateaus. ^^ ^^ The Fermi 
wavelength of several 10 nm for the two-dimensional 
electron gases (see Ref. 24 for instance) is much larger 
than the atomic dimensions. For metals, in contrast, 
the Fermi wavelength is on the order of the interatomic 
distance.^ Therefore, disorder-related effects due to the 
atomic structure are expected to play a more important 
role than for the two-dimensional electron gas devices, 
and material-specific chemical properties enter as an ad- 
ditional aspect. To describe the metallic atomic contacts, 
it is hence necessary to take into account both the atomic 
and the detailed electronic struct ure. l^^r^ ' 

Recent theoretical studies of the thermopower of 
nanocontacts have considered separately the effects of 
disorder and material-specific parameters. Thus, ther- 
moelectric effects for structures with atomically-thin one- 
dimensional wires connected to two-dimensional elec- 
trodes were studied in Refs. 33 34 for various geometries 
and with disorder, but without material-specific parame- 
ters. In contrast in Refs. 135)361 material-specific parame- 
ters were employed, while crucial disorder-related effects, 
resulting in the statistical nature of the experiments,^ 
were neglected by treating only ideal wire geometries. 

Based on molecular dynamics (MD) simulations of 
many junction stretching processes, we study here the 
thermopower of atomic contacts for three different met- 
als, namely Ag, Au, and Pt. We use an atomistic model, 
a material-specific description of the electronic structure. 



and address the influence of disorder at the narrowest 
part of the contacts. The results for the monovalent met- 
als, especially Au, are compared to the available experi- 
mental data,^^ while predictions are made for the multi- 
valent Pt. 

The manuscript is organized as follows. In Sec. [TI|we 
present the methodological details of our atomistic sim- 
ulations and then discuss the results in Sec. IIIII Con- 



clusions and an outlook are given in Sec. |IVl Additional 
material, such as further examples of junction stretch- 
ing processes with animations or technical details, are 
provided in the appendices and in the electronic supple- 
mental material.!^ 



II. THEORETICAL APPROACH 



Our calculations proceed along the lines of Refs.'30'31' 
which we extend here to study the thermopower. The 
calculations can be divided essentially into two parts, 
namely the generation of contact geometries and the de- 
termination of transport properties. 



A. Molecular dynamics simulations 

The contact geometries used in this work for Ag, Au, 
and Pt are those employed in our previous studies. ^^SEil 
For the sake of completeness, we describe their construc- 
tion briefly. 

We determine contact geometries by performing classi- 
cal MD simulations of junction stretching processes. As 
displayed in Figs. [T]to|3) we choose nanowire geometries 
with a central wire (CW) connected to larger-diameter 
electrodes which are attached at the top and bottom. 
Initially, we assume all atoms to be located at the po- 
sitions of a perfect fee lattice, with the (100) direction 
oriented along the z axis. The lattice constants are de- 
termined by minimizing the potential energy of a crys- 
tal. When we use the semiempirical potentials of our MD 
calculations,^^ we obtain lattice constants of 0.408, 0.406, 
0.393 nm for Ag, Au, Pt, respectively, which are close to 
the experimental values."^ The CW consists of 112 atoms 
with 14 layers along the z direction and 8 atoms per layer. 
Its initial length amounts to 2.65, 2.64, 2.55 nm for Ag, 
Au, Pt, respectively, as indicated in Figs.[T]to|3] In each 
case the electrodes at the top and bottom contain 288 
atoms. 

Junction stretching is performed by separating the oth- 
erwise fixed electrodes symmetrically by a constant dis- 
tance during every time step of 1.4 fs. In this process, 
we use periodic boundary conditions along the z direc- 
tion and the minimum image convention for the super- 
cells perpendicular to it .12^ The constant stretching speed 
amounts to 2 m/s. The forces on the wire atoms are cal- 
culated from semiempirical potentials,'^ while their aver- 
age temperature remains at 4.2 K by use of a Nose-Hover 
thermostat. ^^ Different junction evolutions are obtained 



by choosing random starting velocities for the atoms in 
the CW. Every 1.4 ps a contact geometry is recorded, 
and the strain force is calculated.!^ 



B. Transport properties 

Transport properties are determined within the 
Landauer-Biittiker formalism.H^ESESHUmi r^^^^ conduc- 
tance G and thermopower S are expressed as 



G = Goi^o(T), 



K^(T) 



eTK^iry 



(1) 



(2) 



with Kn{T) = ^ dE{E - ^f T{E)[-dEf{E,T)], the 
quantum of conductance Gq = 2e^//i, the absolute value 
of the electron charge e = |e|, the transmission function 
t{E), the Fermi function f{E,T) = {exp[{E-fi)/kBT]^ 
1}~^, the Boltzmann constant /cb, and the electrochem- 
ical potential /i, which approximately equals the metal 
Fermi energy /i ^ Ejr. 

In order to compare to the experiments of Ref. 17, we 
assume in the following a temperature of T = 12 K for 
the determination of the transport properties. While we 
evaluate S via Eq. ([2|, for the conductance the simpler 
low-temperature formula 



G = GoriEp) = Go^TniEp) 



(3) 



yields a good approximation. Here, r has been resolved 
into the contributions of transmission probabilities r^ 
from the individual transmission eigenchannels n. As 
discussed below, they provide important information for 
the understanding of the results. 

We obtain the transmission and the channel decompo- 
sition by use of Green's function techniques, as described 
in Ref. 31, We assume the central device region to be 
equal to the CW, while the remaining atoms of the con- 
tacts are attributed to the electrodes. However, while in 
Ref. |3T] we assumed the unperturbed electrode Green's 
functions to be those of the bulk, in the present study 
we use surface Green's functions. They are determined 
via a decimation technique by modeling the electrodes as 
surfaces of ideal, semi-infinite crystalsJ ^^ * ^^ * Even though 
the results for these two different procedures differ only 
slight ly,S31 surface Green's functions describe the physical 
situation more accurately. 

To perform the energy integrations in Eq. Q, we 
have computed r{E) every 5.6 ps for 11 equally spaced 
points in the energy interval [Ep — A, Ef + A] with A = 
0.05 eV. Kn{T) is then obtained by the integration of 
{E — ii)^ t{E) [—dEf{E^ T)], using a linear interpolation 
for r{E) between the energy sampling points. 

The effective single-particle Hamiltonian and over- 
lap matrices for the evaluation of the transmis- 
sion are obtained from a Slater- Kostei^^ tight-binding 



descriptionpSI^ supplemented by a local charge neutral- 
ity condition.^ Even if this approach is not at the level 
of self-consistent ab-initio methods, it is still atomistic 
and takes into account the symmetries of the 5, p, and 
d valence orbitals for these monoatomic systems. In the 
past, it was used successfully to describe the conduction 
properties of various metallic atomic contacts.^^^^^^^^^J^ 



III. RESULTS AND DISCUSSION 
A. Junction stretching events 

In order to analyze the behavior of the thermopower 
for the metallic atomic contacts, we have simulated 50 
stretching processes for each of the Ag, Au, and Pt 
nanowires. Beside the thermopower we have analyzed 
the strain force and the conductance with its decomposi- 
tion into individual channel contributions. Examples of 
stretching events, which we will discuss in the following 
paragraphs, are shown in Figs. [T] to [3) In all the selected 
cases the contacts break after a dimer contact has formed, 
i.e. a two-atom chain. Further examples of stretching 
events, including the formation of longer atomic chains 
for Au and Pt contacts, are provided in Appendix [A} 

The evolution of the mechanical and transport prop- 
erties for a Ag contact are displayed in Fig. [T] From the 
strain force, shown in the upper panel of the figure, elastic 
and plastic stages can be distinguished. While the force 
increases in a linear manner in the elastic ones, it drops 
suddenly when bonds break and atoms rearrange during 
the short plastic stages. ^^ Similar to what is shown in the 
experimental plots for Au in Ref. dT^ the thermopower 
in the middle panel behaves in a step- wise manner. It 
takes both positive and negative values, and fiuctuates 
around zero. As visible from the dotted lines in Fig. [T] 
the steps in the thermopower typically coincide with plas- 
tic deformations of the contact with smooth changes in 
between. All these features and also the absolute values 
of the thermopower are in agreement with the experi- 
mental results. We note that while data are presented 
in Ref. 17 only for Au, it is stated there that studies of 
Ag and Cu samples showed similar results, which justi- 
fies our comparison. The conductance is displayed in the 
lowest panel of Fig. [T] After an initial increase, it drops 
in a gradual manner. Before contact rupture it decreases 
from a value of around 2Go to a clear plateau at IGq, 
when a single atom is at the narrowest constriction. The 
contact breaks quickly after a dimer has formed. For the 
single-atom and dimer configurations the current is car- 
ried by a single fully t ransmitting channel, as expected 
for the s-valent metals .' ^^ ' ^^ ' ^" ^1^2 

The evolution of junction properties for a Au contact is 
shown in Fig. [2] As suggested by the comparable valence 
electronic structure of the atoms, the results resemble 
those of Ag. The thermopower behaves in a step-wise 
fashion and exhibits gradual changes or plateaus dur- 
ing the elastic stages, interrupted by the plastic defor- 




Figure 1: (Color online) Formation of an Ag dimer contact. 
The upper, middle, and lower panel shows, respectively, the 
tensile force, the thermopower, and the conductance as a func- 
tion of the electrode displacement. For the conductance also 
the decomposition into conduction eigenchannels GoTn is dis- 
played. Horizontal dashed lines indicate the zero level of the 
thermopower or a conductance of ICo, while the vertical dot- 
ted lines mark the main plastic reorganizations of the contact. 
Above and below the graphs snapshots of the contact geom- 
etry during the stretching process are shown, and the length 
of the central wire is given for zero electrode displacement. 



mations. Again it takes both positive and negative val- 
ues, and the thermopower is comparable in magnitude 
to those of Ag and the data presented in Ref. |T7| The 
conductance falls with a revival to a value of I.QGq. It 
decreases then in a rather continuous way to a value of 
1.3Go before it drops to a value slightly below IGq, car- 
ried by a single channel. The conductance changes only 
little, when the single- atom contact transforms into the 
rather stable dimer contact. 

For multivalent Pt the situation is expected to change. 
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Figure 2: (Color online) The same as Fig. IT] but for the for- 
mation of an Au dimer contact. 



and an example is shown in Fig. |3] We observe that also 
here the thermopower fluctuates around zero, but with 
absolute values roughly one order of magnitude larger 
as compared to Ag and Au. The occurrence of these 
strong fluctuations, seen also in the conductance, is due 
to the important contribution of d states at Ef- They 
give rise to a pronounced energy dependence of the trans- 
mission especially for the narrow contacts, as discussed 
in Ref. [SL Also, due to the higher anisotropy of these 
states as compared to the spherically symmetric s va- 
lence orbitals of Ag and Au the junction transport prop- 
erties exhibit a high sensitivity to changes in the atomic 
positions. ^^ While we are not aware of any measurements 
of the thermopower for atomic Pt contacts, we observe 
that the conductance for the single-atom and dimer con- 
tacts is at around 1.5Go (see the electrode displacements 
above 1.02 nm). The multivalent electronic structure and 
hence the contribution of several channels to the trans- 
mission at Ef typically lead to conductance values above 
IGo for one- atom- wide constrictions, in agreement with 
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Figure 3: (Color online) The same as Fig. [T] but for the for- 
mation of a Pt dimer contact. 



the interpretation of measured opening traces and con- 
ductance histograms .I^^IS^H^ 



B. Statistical analysis of atomistic simulations 

Collecting the data of many stretching processes, 
we perform a statistical analysis similar to the 
experiments.^^ For this purpose, we use the 50 simulated 
opening events of the nanowires for each metal. Typi- 
cally the contacts break via the single atom (Ag: 12%; 
Au: 8%; Pt: 4%) or dimer (Ag: 82%; Au: 68%; Pt: 60%) 
configurations mentioned before, with a single atom or a 
chain of two atoms in the narrowest part, respectively. As 
visible from the percentage of occurrence, the dimers are 
generally preferred. However, we find also atomic chain 
geometries with chain lengths of three atoms or more 
(Ag: 6%; Au: 24%; Pt: 36%). For Ag such structures 



are rare. Out of the 50 stretching events, we find at rup- 
ture only three chain configurations with chain lengths 
varying from three to five atoms. While such structures 
are generally believed to be unlikely, ^^ ^^ they have been 
observed in transmission electron microscopy studies.'^ 
For Au and Pt atomic chains occur more frequently in 
our simulations, and they are longer. We have decided to 
exclude the junctions breaking with chains of 10 atoms 
or more in length from our analysis, since experiments 
indicate maximum lengths of up to eight atoms. ^^ ^^ ^^ 
For Au we have additionally excluded an opening pro- 
cess leading to a chain with six atoms, because the ther- 
mopower exhibits peculiar features before contact rup- 
ture related to d states. For this reason we present below 
the statistical analysis of the transport properties with 
50 stretching events for Ag, 46 for Au, and 41 for Pt. 

The thermopower-conductance {S-G) plots are pre- 
sented in Fig. p] The main panels of Fig. pFa-c) show 
all the computed data points. This representation pro- 
vides an impression of the overall scatter. In the insets 
the average value of the thermopower {S)g and the stan- 
dard deviation as = ^J {S'^)g — {S)q are plotted. While 
{S)q is very close to zero for Pt, we find a trend towards 
negative values for Ag and Au. From the theory of the 
free-electron gas in a hyperbolic constriction^^ a ther- 
mopower with a unique, negative sign is expected j ^^ * ^^ * 
We cannot exclude such an origin for the trend, but our 
results with thermopower values of both positive and neg- 
ative sign demonstrate that it is necessary to go beyond 
these simple free-electron-gas models with idealized wire 
geometries to describe even the s-valent metallic atomic 
contacts. Furthermore, within the standard deviation 
the results for {S)q are consistent with a vanishing ther- 
mopower. From Fig.H[a-c) we also observe that the vari- 
ations of S tend to increase for a decreasing minimum 
cross-section of the contacts. Consistent with the sam- 
ple opening traces discussed above, the variations of the 
thermopower around zero are predicted to be one order 
of magnitude larger for Pt than for Ag and Au. 

With respect to {S)q ^ and increasing as for de- 
creasing G, the S-G scatter plot for Au closely resembles 
the experimental one of Ref. [TTl However, the absolute 
magnitude of the scatter of S in our calculations is too 
small by a factor of around 3, and we wil l discuss possible 
reasons for this discrepancy below.^^^^ 

In order to reveal further differences between the 
mono- and multivalent metals, we show the correspond- 
ing S-G density plots in Fig.[4]^d-f). For Ag and Au they 
exhibit a pronounced maximum at (G^S) = (lGo,0), 
while there is no particular feature for Pt at this posi- 
tion. The peak for the 5-valent metals seems plausible. 
In the conductance histograms of Ag and Au pronounced 
maxima occur at IGq-^^ Our calculated conductance 
histograms, shown in Fig.p[a,b), are consistent with this 
experimental finding . ^^ ^^ Since the mean thermopower 
vanishes, maxima at (lGo,0) are thus expected in the 
density plots of Fig. [4Fd,e). The experimental conduc- 
tance histogram for Pt shows instead a rather broad dis- 



tribution between 1 and 2Go with a maximum at around 
1.5Go.^^ ^^ These features are not reproduced in detail 
in our calculations for Pt in Fig. Isjc), which shows a 
broad distribution of frequently occurring conductances 
between 1 and 3Go without a clear maximum. Possible 
reasons for the deviations, such as the limited ensemble 
of considered junction geometries or the approximations 
inherent to our method, have been discussed in Ref. 3TJ 
Still, from the broad distribution of conductance val- 
ues no sharp peak feature at a single conductance value 
is expected in the S-G density plot, in agreement with 
Fig.gf). 

While our discussion of Fig.[4]^d-f) considered until now 
the conductance values, the sharpness of the peak fea- 
ture at (lGo,0) with respect to the thermopower axis is 
somewhat unexpected for Ag and Au. Further insight 
is obtained by plotting the standard deviation as as a 
function of the conductance, as displayed in Fig. |5] The 
calculations show indeed a suppression of as at IGo for 
Ag and Au, which explains the sharpness of the feature in 
the corresponding S-G density plots in Fig. [4]^d,e). Sup- 
pressions of as at higher conductance values are absent 
in Fig. [5|a,b), while shallow minima were measured for 
Au in that region.!^ Consistent with the underestima- 
tion of the scatter of S in Fig. |4je), as for Au is smaller 
than in the experiments by a factor of 3. For the mul- 
tivalent metal Pt no particular features arise through- 
out the whole range of conductance values considered in 
Fig. [sjc), especially there is no minimum of as at IGq. 

In order to describe the behavior of the thermopower 
fluctuations in Fig. [5] in simple terms, we construct an 
"extended single- level model" (ESLM). We assume that 
the transmission is given as 



r(^) = ^r„(E,4"\ri"\r(,"' 



(4) 



with 






p(^)p(^) 

L R 



{E-e^-^y 



p(-) ,p(r.) 

^ T, ^^ 7=? 



(5) 

The function Tn{E^ Cq , F}^ , F)^^) describes a resonance 
in the energy-dependent transmission probability of the 
nth eigenchannel, which is located close to the Fermi en- 
ergy. The resonances arise from the atomic level struc- 
ture at the narrowest part of the contact and quantum 
interference effects related to the disorder in the CW. 
Hence, the "level" Cq specifies the position of the reso- 
nance for the nth eigenchannel, and F)^^ , F)^^ determine 
the width and height of the resonance. For an atomic con- 
tact described by N eigenchannels, we can obtain from 
Eqs. and ([5| the conductance and thermopower us- 
ing the low-temperature expressions of Eqs. ^ and ([2| 
for a given realization of parameters Cq , T^^^ , F]^^ with 
n = 1, . . . , A^. As explained in detail in Appendices[B]and 
ICJ we identify the TniEp^ e^ \ F}^ , F)^ ) with respective 
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Figure 4: (Color online) Scatter plots of the thermopower as a function of the conductance for (a) Ag, (b) Au, and (c) Pt. The 
main panels show all computed data points, while the insets display the average of the thermopower together with its standard 
deviation. Panels (d-f) are the corresponding density plots. For the insets of plots (a-c) and the density plots (d-f), the bin 
size of the conductance is O.lCo- Those of the thermopower amount to 0.075 /iV/K for (d,e) and 0.75 /iV/K for (f). 



eigenchannel transmission probabilities r^ in order to dis- 
tinguish between mono- and multivalent metals. 

As the simplest model for a monovalent metal we as- 
sume that the transmission channels open one by one, 
i.e. that ri = . . . = r^-i = 1 and r^ = G/Gq — (n — 1) 
for n — 1 < G/Gq < n. In this case (see Appendix O) we 
find that 



crs 



OL\/n — G/Gq/u 



(6) 



for G -^ nGo in the conductance interval (n — l)Go ^ 
G < nGo, and a is a constant prefactor of dimension 
V/K. This relation explains the minimum of cr^ at IGq as 
visible in Fig. [5|a,b). Minima at higher conductance val- 
ues nGo with n > 1 are absent for Ag and Au, since the 
assumption that channels open one by one, also know n 
as the "saturation of channel transmission",'^^^^^'^^^ is 
accurately obeyed in our calculations only for the first 
channel, i.e. for conductances below IGq (see Appendix 
[B] and Refs. 30 31). The presence of minima at integer 
conductance values above IGq is hence hard to under- 
stand from our simulations. Frequently occurring junc- 
tion geometries with several atomic chains in parallel^^ 
may provide an explanation, and investigations of junc- 
tions with larger initial diameters could help to resolve 
the puzzle. 

In contrast to the monovalent metals the ESLM of a 
multivalent metal, as presented in Appendix [Cj predicts 
that minima of as are absent, or at least strongly re- 
duced, due to the contribution of several partially open 
transmission channels throughout the whole range of con- 
ductance values. This is consistent with our results from 



the atomistic simulations in Fig. Isjc). 

The behavior of the thermopower has been discussed 
in terms of a different and very detailed model in 
Refs. I17|27t which is applicable to the whole range of 
conductance values. However, apart from the prefactor, 
at G below but close to nGo their expression for as in 
monovalent metals is identical to our Eq. (Ig]). As a com- 
mon feature, both our and their approach explain the 
suppression of the thermopower fluctuations by disorder- 
related quantum interference effects (see also Appendix 
[C| ). In their model the thermopower fluctuations arise 
from the interference of directly transmitted electrons 
and those backscattered elastically in the diffusive re- 
gions in the vicinity of the narrowest part of the contact. 
The backscattering contributions are considered to lowest 
order, and channel transmissions in the ballistic central 
region are assumed to be energy-independent and flxed 
for a certain G. In contrast, in our atomistic simulations 
we determine the energy-dependent transmission of the 
disordered CW, take scattering effects into account to 
all orders, and different realizations of the set of r^ may 
contribute for a given value of the conductance. 

In the experiments both local and nonlocal interference 
contributions will be present and backscattering from de- 
fects up to the coherence length of some 100 nm away 
from the constriction has been reported at the low tem- 
peratures relevant here.!^ Such highly nonlocal interfer- 
ence contributions are clearly not described by our short 
CWs. This may result in an underestimation of the en- 
ergy dependence of the transmission, causing the factor 
of 3 discrepancy with respect to the experimental val- 
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Figure 5: (Color online) Standard deviation as of the ther- 
mopower as a function of the conductance for (a) Ag, (b) 
Au, and (c) Pt. Shown in the background is the conductance 
histogram for each metal, and the number of counts is indi- 
cated on the right y-axis. In each case the bin size of the 
conductance is 0.1 Co- 



nes in Figs. p[b) and p[b) mentioned above.^ Another 
reason for the discrepancy may be the large ensemble of 
experimentally realized contact configurations, not fully 
taken into account in our calculations. 

While measurements of the thermopower and its fluc- 
tuations are not yet available for Pt, Ludoph and van 
Ruitenbeek have shown within their model that as is pro- 
portional to the standard deviation acv of the voltage- 
dependent conductance.^^ The comparison of their re- 
sults for acv for the transition metals niobium and iron 
and the s-valent metals Cu, Ag, and Au shows larger fluc- 
tuations by a factor of two for the former. Our prediction 
of a factor of 10 larger fluctuations for the thermopower 
of Pt as compared with the s-valent metals may be con- 
sidered as an upper bound, and a reduction might arise 
from the finite averaging times in the measurements, for 
instance. 

As another important aspect we have discussed only 
the elastic electronic contribution to the thermopower. 
Inelastic effects due to the electron-phonon coupling may 
lead to modifications, which we expect to be small, how- 
ever, for the experimental conditions considered here. 
Thus, due to its weakness the electron-phonon cou- 
pling can typically be treated perturbatively for metallic 
atomic contacts, ^^ 6Qt65 | 66 ^^^ ^^^ phonon-drag contribu- 



tions should be suppressed because of the small contact 
diameter, the low measure ment te mperatures, and small 
applied thermal gradients j -*^^ * ^^*^^ * Our results, with the 
reasonable agreement between experiment and theory, 
show that the elastic contribution is sufficient to describe 
the main experimental features. 



IV. CONCLUSIONS AND OUTLOOK 

Using molecular dynamics simulations of up to 50 
stretching events, we have analyzed the thermopower 
of atomic contacts of Ag, Au, and Pt. For Ag and 
Au its behavior agrees well with previously reported 
measurements.!^ On a quantitative level, however, the 
experimental scatter of the thern iopow er values is under- 
estimated by a factor of around 3.'^^'^ The thermopower- 
conductance plots show the thermopower to be zero on 
average within the standard deviation for the three met- 
als studied. Furthermore our calculations predict its vari- 
ations around the mean to increase for the narrowest con- 
strictions and to be one order of magnitude larger for the 
multivalent Pt as compared to the s-valent metals Ag and 
Au. At a conductance of one quantum of conductance we 
find, in agreement with the experiment, a suppression of 
the fluctuations of the thermopower for Ag and Au. For 
the multivalent metal Pt possible minima of the fluctu- 
ations of the thermopower should be shifted to higher 
conductance values, but they are predicted to be absent, 
or at least strongly reduced, due to the influence of sev- 
eral partially open transmission eigenchannels. 

Our calculations indicate that the essential characteris- 
tics of the thermopower of metallic atomic-size contacts 
can be understood based on the elastic electronic con- 
tribution combined with effects of rather local disorder 
at the narrowest part of the atomic contacts. Quanti- 
tative differences between our simulations and the ex- 
periment may arise from the larger variability of contact 
configurations in the experiments and distant scatterers, 
not taken into account in the calculations due to the 
limited system size.l^ The quantification of effects due 
to the electron-phonon coupling, which are ex pected to 
be small for metallic single- atom contacts j -*^^ * ^^ "^ consti- 
tutes a challenging task for the future. 
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Figure 6: (Color online) The same as Fig. [T] for another Ag 
contact forming a dimer before rupture. 
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Figure 7: (Color online) The same as Fig. [T] but for a Au 
contact forming a chain of four atoms before rupture. 



Appendix A: Further junction stretching events 

Here, we provide further examples for the evolution of 
junction properties upon stretching for each of the three 
metals studied. For Au and Pt we have selected junctions 
forming chains of several atoms in length before rupture. 
Animations showing the stretching processes displayed in 
Figs.[T]to[3]and Figs. IB] to [8| can be found as supplemental 
material on the web.^^ 

In Fig. [6] we show an Ag contact. After the initial 
stage, the conductance evolves from a value of around 
4Go via various plateaus to around 2Go until a rather 
stable dimer is formed. For the dimer contact the con- 
ductance is pinned closely to IGq and, at the same time, 
the thermopower and its fluctuations are suppressed to 
zero. This particular example agrees well with the ob- 
servation of the small {S)q and minimum of as for Ag, 
as shown in Figs. [4] and [5J when the conductance arises 
from a single completely open conduction channel. 



For Au we display in Fig. [7] a contact forming a four- 
atom chain before breaking. The chain is increasing suc- 
cessively in length, starting from a single atom for elec- 
trode displacements between 0.84 and 1.12 nm. Two of 
them are present from 1.12 to 1.2 nm, three from 1.2 to 
1.51 nm, and four from 1.51 to 1.59 nm. As soon as the 
single-atom contact has formed, only a single channel 
contributes to the conductance. For the one- to three- 
atom chain configurations the conductance is very close 
to IGo and the thermopower varies around zero with a 
relatively small amplitude. When the four-atom chain 
has formed, the conductance is somewhat suppressed and 
the thermopower increases in absolute magnitude. 

In Fig. [8] we show a Pt contact forming a chain of 
six atoms before rupture. As in Fig. [3] the conductance 
and thermopower during the stretching show much larger 
variations than for Ag and Au, and also the typical values 
of the thermopower are around one order of magnitude 
larger than for the 5-valent metals. After the single- atom 
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Figure 8: (Color online) The same as Fig. [T] but for a Pt 
contact forming a chain of six atoms before rupture. 



contact has formed at around 0.66 nm, more atoms are 
incorporated into the chain at the elongations marked 
with dotted lines (i.e. at 0.95, 1.11, 1.36, 1.38, and 1.59 
nm), until the final length of six atoms is reached. In this 
atomic-chain regime we generally observe two to three 
partially open transmission channels, often leading to 
conductances exceeding IGq. In addition, for these nar- 
rowest contacts the thermopower shows its largest fluc- 
tuations with the average remaining close to zero. 



Appendix B: Average eigenchannel transmission 
probabilities 

For the monovalent metals transmission eigenchannels 
are typically assumed to open one by one, while for mul- 
tivalent metals there are several channels contributing 
already at the lowest conductances l^Il^^MJ \Ye can in- 
spect the validity of these assertions by examining the 



average eigenchannel transmission probabilities obtained 
from our MD transp ort sim ulations as a function of the 
junction conduct ance.'^^^ 

Using the results of the present study, we display the 
average eigenchannel transmissions in Fig. |9] We do in- 
deed observe that for the monovalent metals there is es- 
sentially only a single open channel for < G < IGq. 
However, the first channel is not fully open at IGo, i-e. 
does not reach a transmission of one. For G > IGq the 
saturation of channel transmission still seems to be a rea- 
sonable approximation for Au. However, the results for 
Ag rather suggest an interpretation in terms of a simul- 
taneous opening of two 7r-like channels, whose degener- 
acy is lifted in the disordered junction geometries. For 
Pt there are significant contributions from at least two 
partially open channels throughout the whole range of 
conductance values. 



Appendix C: Fluctuations of the thermopoAver 
within the extended single-level model 

1. Extended single- level model 

We consider a model for the charge transport through 
metallic atomic contacts in order to explain the results 
of our atomistic simulations in simple terms. We assume 
that the energy-dependent transmission r{E) arises from 
those of independent single- level transmissions, as de- 
scribed by Eqs. Q and ([5|. Setting the Fermi energy 
to zero, Ep = 0, we can determine the conductance and 
thermopower in the limit of low temperatures as 



G = Got{E)\^^^^ 



GoE 



p(^)p(^) 

L R 



(Cl) 



^(4"0 +(ri") + r(,")) /4 



-So 



-So 



dET{E) 



t{E) 



E=Ef 



2^(n)p(n)p(n) 
L R 






(e^-))V(ri-)+r^-))'/4 



>/G 



(C2) 



with Go = 2eV/i and ^o = ir'^klT/Se. 

During the stretching of a contact the atomic positions 
change, causing related variations in the electronic struc- 
ture and in the quantum interference pattern. Therefore, 
we consider the parameters e^^ , T]^^ , F)^^ to be indepen- 
dent random variables distributed with probability den- 
sities Pe(4''^), Pri^L^). Pri^R^). respectively. For the 
discussion of the fluctuations of the thermopower, we de- 
termine the standard deviation 



as{G) = 7(52)^ 



{Sfo, 



(C3) 
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Figure 9: (Color online) Average eigenchannel transmission probabilities {Tn)Q as a function of the conductance for (a) Ag, 
(b) Au, and (c) Pt, as determined from the atomistic simulations presented in the manuscript. Error bars show the standard 
deviation and the dashed horizontal lines indicate unit average transmission. See also Refs. [30l31l for further explanations. 



where ( )q means a conditional average over all con- 
figurations yielding the conductance G. For a sys- 
tem of N levels the configurations are labeled by x = 



.(1) -r(i) -r(i) 



V^O 5-'- L ' -■- i? ' • • • ' ^0 

probability density 



{N) j.{N) ^(AT) 



, F]^ 7 T)^ ^). By introducing the 



N 



p{x) 



n^.(4"^)^r(ri"))Pr(r(,"^ 



we can express such an average as 

(5')^= f dxS{x)p{x\G) 



(C4) 



(C5) 



with the conditional probability density p{x\G) = 5{G — 
G{x))p{x)/ fdxS{G - G{x))p{x). Here, G{x) and S{x) 
are determined from Eqs. ( |C1[ ) and ( |C2[ ), respectively. 

In order to describe material-specific properties and 
to understand the behavior of as in Fig. [5J we need to 
know how the individual resonances rn{E^ e^^ , T)^^ , F)^^ ) 
of the eigenchannels in Eq. ^ contribute to the trans- 
mission for a given conductance. For this purpose, 
we identify rn(£^F, Cq , F}^^ , F)^^) with the measurable 
eigenchannel transmission probability r^ to complete the 
ESLM. This identification procedure will be explained 
further in the next subsections. 



2. Monovalent metals 

As visible from Fig. [9J the saturation of channel 
transmission is well fulfilled in our simulations onl y for 
G < IGq. However, to compare to the literature^^l^ 
and since we are mainly interested in conductance val- 
ues around IGq, we adopt the following ideal model to 
describe the monovalent metals. For a given G with 
n — 1 < G/Gq < n, we determine the transmission as 



TiE) 



(n) p(n) p(n) 



i + T„(i?,er,rr,r^"0- (C6) 



Hence, n = . . 
1 < G/Go 



n 



< 



= Tn-l = 1, Tn = G/Gq - (n - 1) for 

n. This behavior of the transmission 



probabilities r^ = Tn{EF^ e^ \ T^T ' ^^r!) ^^ ^ function of 
the conductance is visualized in Fig. IToFa). 

In Fig. [5] it is seen that (Js is suppressed at G ?^ IGq 
for Ag and Au. In order to analyze this, we consider 
a conductance G close to but smaller than nGo in the 
ESLM. For (n - 1) < G/Gq < n we obtain 



G/Go 



p(^)p(^) 

L R 



(4«)) +(ri") + r(,")) /4 



S/So 



2e, 



(") 



(G/Go + l-ny 



p(n)p(n) 
L R 



G/G^ 



Since r^ ~ 1 for G ^ ^Gq, Eq. (C7) requires F 



<n) 



Sn) 



(C7) 



(C8) 



<n) 



F)^ and e^ ~ 0. Therefore, we consider the symmetric 

junction model with F = T^ = F)^ > and eo = Cq . 
Since we are dealing with atomic contacts, we can expect 
Pr(r) to be centered at a nonvanishing, positive value, 
while we assume Pe(eo) to be symmetric and maximal at 
eo = 0. With these assumptions, it follows that 



a| = (5%=< 



\ (G/Gq + l-nf 



rVc 



(G/Go)^ 



(n-G/Go), 

(C9) 

where (l/F^)^ approaches a fixed value in the limit G -^ 
uGq. From the expression, we hence obtain Eq. (|6|. 



3. Multivalent metals 

Based on the results in Fig. [9Fc), we want to construct 
an ESLM for the multivalent Pt. To simplify the situ- 
ation, we consider only two partially open channels and 
write 
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Figure 10: (Color online) Behavior of the eigenchannel transmission probabilities Tn and the standard deviation of the ther- 
mopower as as a function of the conductance for (a,c) a monovalent metal and (b,d) a multivalent metal. For (c) and (d) 
we used 7 = 0.6 eV and 0.1 eV, respectively, and the bin size was set to AC = 0.04^0 in the main panels. In the insets of 
panel (c) we have chosen a smaller bin size of AG = 0.002^0 to resolve the selected region more clearly and show the function 
Gs — OL^n — G/Go/n for G < nGo and n = 1, 2 as a red line, assuming a — 1.44 /iV/K. 
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if < a/Go < l/ai 

if 1/ai < G/Go < (2ai - as) /aj ' 



(CIO) 



with n = aiG/Go, ts = asG/Go for < G/Gq < 1/ai, 
and Ti = 1, r2 = a2/ai + ai(G/Go - l/^i), ts = 
a2{G/Go - l/ai) for 1/ai < G/Gq < (2ai-a2)/a{. 
The model requires ai + a2 = 1, and ai = 0.7 seems 
to be a reasonable choi ce for Pt. The upper bound for 
the conductance in Eq. (CIO) considers in each case the 



full opening of a channel, and the model can easily be 
extended to include further partially open channels or 
describe larger values of G. The behavior of the eigen- 
channel transmission probabilities as a function of the 
conductance is visualized in Fig. floFb). 



We note that in the conductance range < G/Gq < 
1/ai, for a given realization Cq ,r]^ ,r]^^ yielding ri = 

Ti{EF,e^o\^L\^R^)^ we guess values for e^'\rf ,rg^ 
until (1 — p)a2Ti/ai < r2 < (1 +p)a2ri/ai with T2 = 

(2) (2) (2) 

T2{Ef^ Cq , F)^ , F)^^) and with a small tolerance param- 
eter p = 0.05. For conductance values G/Gq > 1/ai we 
proceed similarly. 



Statistical analysis of the fluctuations of the 
thermopower 



By treating the parameters Cq , F 



dent random numbers with 



(n) p(n) 
L 

< 



in) 



, F]^^ as indepen- 
0.1 eV and < 



F^"^,F^'^ < 7, we det ermin e (G, S) pairs from the pre- 



script ions in Eqs. ( |C6[ ) and ( |C10| ). We choose a bin size 
AG and perform the typical statistical analysis. Specif- 
ically, we use {S)q = ^jSj/M where Sj are all those 
thermopower values of the pairs with G — AG/2 < Gj < 
G + AG/2, and M is the total number of such pairs. 
Analogously, we obtain the standard deviation of the 

thermopower as crs = xj^ZA^j ~ {^)g)'^ 1^- 

By assuming Cq to vary symmetrically around zero, 
we obtain {S)g ~ in the ESLM. In this way we can- 
not describe the systematic deviations from zero with a 
unique sign, predicted and measured for quantum point 
contacts realized in a two-dimensional electron gas.^^^"^ 
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However, our assumption is consistent with the experi- 
ments on atomic contacts, where (6^)^ was found to scat- 
ter around zero largely within the cr^P^^ 

From the numerical analysis we obtain the (Js{G) 
curves shown in Fig. 10^c,d). In order to obtain stan- 
dard deviations comparable to the atomistic simulations 
we have set the maximum of the couplings to 7 = 0.6 
eV for the monovalent metal. Because of the sharper res- 
onances due to d states for Pt (compare the Figs. 3, 7, 
and 9 in Ref.|3T]) a smaller value 7 = 0.1 eV is neede d for 
the model of the multivalent metal. Based on Eq. ( C9 ) 



we observe that a smaller 7 should generally increase as- 
To obtain ESLM results for as of the same size as in the 
experiments on Au,^*" 7 should correspond to around 0.2 
eV. Since we would expect line widths F^^ , F^^ around 
1 eV for atomic levels, the small values of even the max- 
imum couplings are consistent with the interpretation 
that the resonances rn(£^, Cq , F)^ , F)^^) are not simply 
due to atomic orbitals in the narrowest part of the con- 
tact, but that they arise essentially from disorder-related 
quantum interference effects. 

A clear suppression of the fluctuations of the ther- 
mopower is visible in the regions G ^ nOo in Fig. [lOJ^c) 
for the model of a monovalent metal. The insets of that 
panel illustrate that as = ay^n — G/Go/n at G ^ nGo, in 



as expected from Eq. (|6|. The suppression at IGq is 
consistent with the results of our atomistic simulations 
and the experiments.^ For Ag and Au the dip at IGq 
in Fig. |5|a,b) is reduced in depth as a result of the in- 
complete opening of the first conductance channel, the 
small contribution of further partially open conduction 
channels [see Fig. |9Fa,b)], and the limited statistics with 
the related large AG = O.IGq. These effects are en- 
hanced at the higher conductance values, where no sup- 
pression of cr^ at 2Go is visible in Fig.[5|a,b). Also in the 
experiment!^ only shallow minima have been observed at 
the positions nGo for n > 2. 



For the ESLM of a multivalent metal shown in 
Fig. floFd), the first transmission channel is fully open 
only at G = Go/ai ~ 1.43Go, so that a possible suppres- 
sion of the fluctuations of S is shifted to conductance val- 
ues above IGq. However, the contribution of the second 
partially open channel in the ESLM, and possibly further 
channels in the experiments [see also Fig. |9jc)], strongly 
washes out the expected minimum at Go/ai- Together 
with the incomplete opening of the dominant conduction 
channel [see Fig. [OJ^c)] these effects explain the absence 
of any clear minima of as for the atomistic simulations 

Fig. He). 
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